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TITLE OF THE INVENTION 

METHOD AND SYSTEM FOR VASCULAR 
ELASTOGRAPHY 

5 

FIELD OF THE INVENTION 

The present invention relates to vascular tissue 
characterization. More specifically, the present invention is concerned with 
10 a method and system for vascular elastography imaging. 

BACKGROUND OF THE INVENTION 

In the early nineties, Ophir et a\., (1991) introduced 
15 elastography, which is defined as biological tissue elasticity imaging. 
Primary objectives of elastography were to complement B-mode 
ultrasound as a screening method to detect hard areas in the breast 
[Garra et a/., 1997]. 

20 Within the last few years, elastography has also found 

application in vessel wall characterization using endovascular catheters 
[Brusseau etal., (2001); de Korte et a/., (1997-2000b)]. Indeed, changes 
in vessel wall elasticity may be indicative of vessel pathologies. It is 
known, for example, that the presence of plaque stiffens the vascular wall, 

25 and that the heterogeneity of its composition may lead to plaque rupture 
and thrombosis. As indicated below, the mechanical properties of plaques 
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were also recently studied non-invasively. 

A) Non-Invasive Vascular Elastography (NIVE) 

5 Recently, several groups have proposed different 

approaches to non-invasively characterize superficial arteries by using 
standard extra-corporal array transducers. Namely, Bang et al. (2003) 
developed a method to analyze pulsatile motion of carotid artery plaque 
from a sequence of radio-frequency (RF) images. Kanai et al. (2003) 
10 proposed to compute a map of elastic moduli to characterize the carotid 
artery. Mai and Insana (2002) proposed to monitor deformations in 
tissues surrounding superficial arteries; a tissue-like gelatin elasticity-flow 
phantom and in vivo scanning of the normal brachial artery were used to 
validate the method. 

15 

However, most of the conventional methods in 
elastography only provide the map of the strain distribution in the direction 
of the ultrasound beam propagation (axial strain, or radial strain in 
endovascular elastography). Tissue stiffness is conventionally 
20 represented with a color code where dark and bright are associated to 
hard and soft tissues. This can set a potential limitation in non-invasive 
characterization of vessel walls with an extracorporal ultrasound probe, 
since the ultrasound beam propagates axially whereas the tissue motion 
runs radially. 



25 
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Indeed, it has been shown that motion parameters might 
be difficult to interpret when tissue motion and the ultrasound beam do 
not occur in the same orientation. As a consequence of that, the axial and 
lateral strain parameters are subjected to hardening and softening 
5 artifacts, which are to be counteracted to appropriately characterize the 
vessel wall. These artifacts are also known as mechanical artifacts that 
can be defined as a wrong association of a strain pattern to hard or soft 
tissues. Mechanical artifacts originate from a combination of the intrinsic 
mechanical properties of the underlying tissue, its geometry, and its 
10 kinetics. 

B) Non-Invasive Micro-Vascular Elastography (MicroNive) 

The following literature review is focusing on the effect of 
15 hypertension on the remodeling of the vascular wall. However, the 
proposed technology is not restricted to this application and concerns 
imaging of the mechanical structures of small vessels in humans and 
small animals such as rats and mice. The targeted diseases are not 
restricted to hypertension and include any pathology affecting the 
20 mechanical properties and structures of the vascular wall such as 
atherosclerosis, for which specific animal models were developed. 

Regarding the investigation of the phenotyping in 
hypertention (HT) with genetically-engineered rat models, it has been 
25 shown that the structural and mechanical properties of the arterial tree are 
altered. Many researchers have investigated such an assertion. For 
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example, Intengan et al. (1998a, 1998b), using DOCA-salt rats, 
specifically studied the interaction of vasopressin and endothelin 1 (ET-1 ) 
in the pathogenesis of the structural vascular alterations. Their results 
demonstrated that the DOCA-salt model of HT is associated with vascular 
5 growth (increased media width, media-lumen ratio, and a growth index of 
44%) in the absence of changes in vascular distensibility. On the other 
hand, many investigations have also been conducted on the role of the 
proximal arteries in HT. Namely; Cantini et al. (2001 ) observed an aortic 
wall stiffness increasing with age in Wistar rats. Tatchum-Talom et al. 

10 (2001) demonstrated that the aortic stiffness is decreased in estrogen- 
deficient rats. Many other researches demonstrated the implication of the 
proximal arteries in the pathophysiology of HT [Johns et al., (1998); Si et 
al., (1999); Goud ef a/., (1998); Zhao ef a/., (2002)]. However, those 
experiments were conducted ex vivo and required the sacrifice of the 

15 animals. Since the most relevant insights into vascular diseases should 
come from in situ investigations, there is a need for non-invasive micro- 
vascular ultrasound elastography (MicroNIVE). The availability of such a 
mechanical characterization non-invasively could lead to significant new 
discoveries in functional genomics and pharmacogenetics [Hamet et al., 

20 (2002)]. 

C) Endovascular Elastography (EVE) 

Atherosclerosis, which is a disease of the intima layer of 
25 arteries, remains a major cause of mortality in western countries. This 
pathology is characterized by a focal accumulation of lipids, complex 



WO 2005/074804 



PCT/CA2005/000162 



5 



carbohydrates, blood cells, fibrous tissues and calcified deposits, forming 
a plaque that thickens and hardens the arterial wall. A severe 
complication of atherosclerosis is thrombosis, a consequence to plaque 
rupture or fissure, which might lead, according to the event localization, to 
5 unstable angina, brain or myocardial infarction, and sudden ischemic 
death [Falk, (1989); Davies and Thomas (1985); Zaman et a/., (2000)]. 
Plaque rupture is a complicated mechanical process, correlated with 
plaque morphology, composition, mechanical properties and with the 
blood pressure and its long term repetitive cycle [Fung, (1993); Falk, 
1 0 (1 992)]. Extracting information on the plaque local mechanical properties 
and on the surrounding tissues may thus reveal relevant features about 
plaque vulnerability [Fisher et a/., (2000); Ohayon et a/., (2001)]. 
Unfortunately no imaging modality, currently in clinical use, allows the 
access to these properties. 

15 

So far, diagnosis and prognosis of atherosclerosis 
evolution in humans mainly rely on plaque morphology and vessel 
stenosis degree. This information can accurately be accessed with 
Intravascular UltraSound (IVUS) imaging, since this modality provides 

20 high resolution cross-sectional images of arteries. Accurate quantitative 
analysis of the disease is thus easily performed by precise measurements 
of the lumen area, arterial dimensions and dimensions specific to the 
plaque. Moreover, IVUS permits the qualitative characterization of plaque 
components, but roughly, in terms of fatty, fibrous or calcified plaques and 

25 with possible misinterpretations. This makes IVUS, alone, insufficient to 
predict the plaque mechanical behavior. However, elastic properties of 
vessel walls can be derived from radio-frequency (RF) or alternatively 
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from B-mode IVUS images, by integrating elastographic processing 
methods. Indeed, endovascular ultrasound elastography (EVE) is an in- 
development imaging technique that aims to outline elastic properties of 
vessel walls. Its principle consists of acquiring sequences of cross- 
5 sectional vessel ultrasound images, while the vascular tissue is 
compressed by applying a force from within the lumen. Strain distribution 
is then estimated by tracking, from the signals, the modifications induced 
by the stress application. In practice, in EVE, such a stress can be 
induced by the normal cardiac pulsation or by using a compliant 
1 0 intravascular angioplasty balloon. 

Several approaches have been proposed to assess 
tissue motion in EVE. Whereas one-dimensional (1D) motion estimators 
are likely more sensitive to pre- and post-motion signal decoherence, two- 

15 dimensional (2D) motion estimators are expected to be more reliable. 
However, the most commonly used motion estimators in EVE applications 
are 1 D correlation-based techniques. This choice is mainly dictated by the 
ability of such estimators to be implemented; they also may provide real- 
time tissue motion estimates. In 1D correlation-based tissue motion 

20 estimators, the displacement between pre- and post-motion pairs of RF- 
or B-mode lines is determined using cross-correlation analysis. This 
technique was used to investigate EVE feasibility on vessel-mimicking 
phantoms [de Korte et a/., (1997)], on excised human femoral and 
coronary arteries [de Korte etal., (1998); (2000a)], and in vivo on human 

25 coronary arteries [de Korte et al. % (2000b)]. 
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A different implementation of the 1D cross-correlation 
technique was proposed by Brusseau et al. (2001 ) to investigate a post- 
mortem human excised carotid artery. Brusseau computed an adaptive 
and iterative estimation of local scaling factors, using the phase 
5 information between pre- and post-compression RF signals. The authors 
suggested that this approach may be less sensitive to decorrelation noise 
than conventional 1D correlation-based estimators. On the other hand, 
others also proposed to assess local scaling factors, but in the frequency 
domain [Thalami et a/., (1994)]. They presented some initial in vitro and in 
10 vivo results that were obtained with this spectral tissue strain estimator. 
Envelope B-mode data were used in this last study. However, no further 
validation of the spectral approach was so far conducted in EVE. 



In vivo applications of EVE are subjected to many 
15 difficulties. For instance, the position of the catheter in the lumen is 
generally neither in the center nor parallel to the vessel axis, and the 
lumen geometry is generally not circular. In such conditions, tissue 
displacements may be misaligned with the ultrasound beam, introducing 
substantial decorrelation between the pre- and the post-tissue- 
20 compression signals. In addition, although the ultrasound beam 
propagates close to parallel with the tissue motion in EVE, providing the 
full strain tensor should improve the characterization of complex 
heterogeneous tissue structures that may deform unpredictably following 
the cardiac pulsation of the vessel. The complex heterogeneous nature of 
25 plaques may indeed induce 1D decorrelation due to the complex 3D 
movement of the tissue structures. Regarding that, 1 D estimators may not 
be optimal if such decorrelation is not appropriately compensated for. 
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Ryan and Foster (1997) then proposed to use a 2D correlation-based 
speckle tracking method to compute vascular elastograms. This approach 
was experimented on envelope B-mode data from in vitro vessel- 
mimicking phantoms. No further validation was however conducted by this 
5 group. 



Another potential difficulty, that is associated with EVE in 
vivo applications, stems from the eventual cyclic catheter movement in 
the vessel lumen. Owing to the pulsatile blood flow motion, catheter 

1 0 instability may constitute another source of signal decorrelation between 
pre- and post-compression signals. To that, Shapo etal. (1996a; 1996b) 
proposed the use of an angioplasty balloon to stabilize the catheter in the 
vessel lumen. Tissue motion was assessed using a 2D correlation-based 
phase sensitive speckle tracking technique. Preliminary results from 

15 simulations and from in vitro vessel-mimicking phantom investigations 
were presented; envelope B-mode data were used. 



OBJECTS OF THE INVENTION 



20 An object of the present invention is therefore to provide an 

improved method and system for vascular elastography. Another object is 
to provide a method and system to non-invasively map the elastic 
properties of vessels. 
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SUMMARY OF THE INVENTION 

5 According to a first aspect of the present invention, there is 

provided a method for vascular elastography comprising: 

providing pre-tissue-motion and post-tissue-motion images in 
digital form of a vessel delimited by a vascular wall; the pre-tissue-motion 
and post-tissue-motion images being representative of first and second 
10 time-delayed configuration of the vessel; 

partitioning at least portions of both the pre-tissue-motion 
and post-tissue-motion images into corresponding data windows; 

approximating a trajectory between the pre-tissue-motion 
and post-tissue-motion images for corresponding data windows; and 
1 5 using the trajectory for each data window to compute a strain 

tensor in each data window. 

The method can be adapted for non-invasive vascular 
ultrasound elastography (NIVE) to non-invasively characterize superficial 
20 vessels such as carotid, femoral arteries, etc. NIVE is of clinical values for 
the purpose of diagnosis and follow-up of vascular pathologies. 

The method can further be adapted for non-invasive vascular 
ultrasound micro-elastography (MicroNIVE) for characterizing small 
25 superficial vessels in humans and animals. More specifically but not 
exclusively, MicroNIVE is of value in functional genomics to investigate 



WO 2005/074804 



PCT/CA2005/000162 



10 

phenotyping in hypertension with genetically-engineered rat models. 

The method for vascular elastography according to the first 
aspect of the present invention can also be adapted for endovascular 
5 ultrasound elastography (EVE) for invasive characterization of vessels 
using catheter-based techniques. More specifically but not exclusively, 
EVE is used to investigate coronary diseases in humans. 

The method for vascular elastography according to the first 
10 aspect of the present invention can also be adapted to other imaging 
technologies such as, but not exclusively, to magnetic resonance imaging 
(MRI), optical coherence tomography (OCT) or Doppler-based ultrasound 
imaging for the non-invasive and invasive characterization of vessels, 
providing that the imaging techniques can provide the assessment of 
15 tissue motion. 

According to a second aspect of the present invention, there 
is provided a system for vascular elastography comprising: 

an ultrasound system for acquiring pre-tissue motion and 
20 post-tissue motion radio-frequency (RF) images of a vessel; the pre-tissue 
motion and post-tissue motion images being representative of first and 
second time-delayed configuration of the vessel; 

a controller, coupled to the ultrasound system, i) for receiving 
the pre-tissue motion and post-tissue motion RF images, ii) for digitizing 
25 the pre-tissue motion and post-tissue motion RF images, iii) for partitioning 
both the pre-tissue motion and post-tissue motion RF images within the 
vascular wall into corresponding data windows, iv) for approximating a 
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trajectory for each the data windows; and v) for using the trajectory for 
each the data window to compute a strain tensor in each data window; 
and 

an output device coupled to the controller to output 
5 information related to the strain tensor in each data window. 

Other objects, advantages and features of the present 
invention will become more apparent upon reading the following non 
restrictive description of preferred embodiments thereof, given by way of 

10 example only with reference to the accompanying drawings. It is to be 
noted that the examples presented hereafter were based on the analysis 
of the ultrasound RF signals. The present invention is not restricted to RF 
or B-mode signals and may be applied to any new ultrasound modalities 
providing tissue movements. The RF signals may be seen as the raw data 

15 from which all current imaging modalities available on the market were 
developed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 In the appended drawings: 

Figure 1 is a block diagram of a system for vascular 
elastography according to a first illustrative embodiment of a first aspect of 
the present invention; 

25 

Figures 2 and 3 are respectively a flowchart and a block 
diagram illustrating a method for vascular elastography according to a first 
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illustrative embodiment of a second aspect of the present invention; 

Figure 4 is a schematic view illustrating a two-dimensional 
partitioning of a region of interest (ROI) within a vascular wall, part of the 
5 method illustrated in Figures 2 and 3; 

Figure 5 is a block diagram illustrating a method for vascular 
elastography according to a second illustrative embodiment of the first 
aspect of the present invention; 

10 

Figures 6A-6F are theoretical gray-scaled displacement 
fields and elastograms illustrating motion parameters for a pressurized 
thick-wall cylindrical blood vessel, embedded in an elastic infinite medium; 

15 Figures 7A-7E are theoretical gray-scaled displacement 

fields and elastograms illustrating radial strain and strain decay for a 
homogeneous vessel wall; 

Figures 8A-8C are respectively gray-scaled elastograms (8A- 
20 8B) obtained and a graph illustrating the comparison between the radial 
strain from Figures 7 and the Von Mises (VM) parameter ; 

Figure 9 is a schematic view of an experimental set-up used 
to produce mechanical deformation of polyvinyl alcohol cryogel (PVA-C) 
25 vessel-mimicking phantoms, and to collect RF ultrasound data 
incorporating the system from Figure 1; 
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Figure 1 0 is a schematic view of the vascular flow phantom 
from the experimental set-up from Figure 9; 

Figures 11A-11C are schematic views of the moulds that 
5 were used to construct the double-layer PVA-C vessel from Figure 10; 

Figures 12A-12C are respectively a B-mode image, a Von 
Mises (VM or £,) elastogram obtained using the method from Figure 2 and 
the set-up from Figure 9 and a graph illustrating the average of 5 axial 
10 lines chosen in the middle of E, in the Figure 12B; Figure 12A being labeled 
"Prior Art"; 

Figures 13A-13B, which are labeled "Prior art", are 
respectively B-mode image of a carotid artery acquired from a healthy 
15 volunteer, and a manually segmented B-mode image of the vessel wall; 

Figures 13C-13D, are gray-scaled elastograms computed 
from data acquired at two different locations of the carotid artery from 
Figures 13A-13B, using the method from Figure 2; 

20 

Figures 14A-14B, which are labelled "Prior Art", are B-mode 
images acquired over longitudinal sections of the carotid artery of 
respectively a normotensive and a hypertensive rat; 

25 Figures 14C-14H are axial strain "gray-scaled" elastograms 

of the carotid artery of six different rats, three normotensive (C-E) and 
three hypertensive (F-H) obtained using the method from Figure 2; 
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Figure 15 is a schematic view illustrating the image 
acquisition process part of a method for endovascular elastography 
according to a third illustrative embodiment of the second aspect of the 
5 present invention; 

Figure 1 6 is a schematic view illustrating an "ideal" plaque in 
a vascular tissue representation; 

10 Figure 17A-17B are respectively an in vivo intravascular 

ultrasound cross-sectional image of a coronary plaque and a two- 
dimensional finite element mesh of the unload real geometry with spatial 
distribution of the constituents from the plaque from Figure 17A; Figure 
17A being labeled "Prior Art"; 

15 

Figures 18A-18D are respectively a theoretical "gray-scaled" 
elastogram of a radial strain computed for an idealized plaque; a graph 
illustrating theoretical radial strain distributions taken along the respective 
lines from Figure 18A; a radial strain "gray-scaled" elastogram obtained 
20 using the endovascular elastography method according to the third 
illustrative embodiment of the second aspect of the present invention; and 
a graph illustrating the radial strain distributions taken along the respective 
lines from Figure 18C; 

25 Figures 19A-19C are respectively a strain-decay- 

compensated "gray-scaled" elastogram obtained using the endovascular 
elastography method according to the third illustrative embodiment of the 
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second aspect of the present invention; and one-dimensional vertical and 
horizontal graphs taken along the respective lines from Figure 19A; 

Figures 20A-20C are respectively a theoretical radial strain 
5 elastogram of the coronary artery illustrated in Figure 17A; and one- 
dimensional vertical and horizontal graphs taken along the respective lines 
from Figure 20A; 

Figures 21A-21C are respectively radial strain "gray-scaled" 
elastogram computed for the coronary artery illustrated in Figure 17A 
using the method for endovascular elastography according to the third 
illustrative embodiment of the second aspect of the present invention; and 
one-dimensional vertical and horizontal graphs taken along the respective 
lines from Figure 21 A; 

Figures 22A-22C are respectively a strain-decay- 
compensated "gray-scaled" elastogram of the coronary artery illustrated in 
Figure 17A obtained using the endovascular elastography method 
according to a third illustrative embodiment of the second aspect of the 
present invention; and one-dimensional vertical and horizontal graphs 
taken along the respective lines from Figure 22A; 

Figure 23 is a schematic view of an experimental set-up 
including a system for endovascular elastography according to a second 
25 embodiment of the first aspect of the present invention; 



10 



15 



20 



Figures 24A-24C, which are labelled "Prior Art", are 
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respectively a histological section of a post-mortem excised human carotid 
artery with a very thin plaque; a close-up view of the atherosclerotic region 
taken from Figure 24A; and a log-compressed IVUS image of the carotid 
section; and 

5 

Figures 25A-25J are "gray-scaled" elastograms computed for 
consecutive increasing physiologic fluid pressure levels for the carotid 
artery illustrated in Figures 24A-24C using the method for endovascular 
elastography according to the third illustrative embodiment of the present 
10 invention. 

DETAILED DESCRIPTION OF THE INVENTION 

15 A system 10 for vascular elastography according to a first 

embodiment of a first aspect of the present invention will now be 
described with reference to Figure 1. More specifically, the system 10 
allows for non-invasively characterizing arteries. Whereas not restricted to, 
this system allows predicting risks of vascular tissue rupture due to the 

20 presence of atherosclerotic plaques and potentially vascular aneurysms. 
Since vascular tissue rupture due to atherosclerotic plaques and 
aneurysms is believed to be well known in the art, it will not be described 
herein in more detail. 

25 The system 1 0 comprises an ultrasound system 1 1 including 

an ultrasound instrument 12 provided with a scanhead 20 including an 
ultrasound transducer. The instrument 12 is coupled to an analog-to-digital 
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acquisition board 14 of a controller 16 via a radio-frequency (RF) pre- 
amplifier 18. 

For NIVE, the ultrasound instrument 12 is configured for 
5 extracorporal measurement, while for MicroNive, it is in the form of an 
ultrasound biomicroscope. The ultrasound system 1 1 is configured with 
access to RF data so as to allow computing vascular elastograms of 
vessels. Examples of such ultrasound system 1 1 are the ES500RP from 
Ultrasonix for NIVE, and the high-resolution VS-40 or Vevo660 from 
10 Visualsonics for MicroNive. An ultrasound system from another type or 
having other configurations can also be used. 

The ultrasound instrument 12 provides an RF output from 
which the received RF data were transferred to the pre-amplifier 18. An 
1 5 example of pre-amplifier that can be used is the Panametrics, model 5900 
PR. Of course, other pre-amplifier can alternatively be used. 

The acquisition board 14 allows digitizing the pre-amplified 
signals from the pre-amplifier 18. An example of acquisition board is the 
20 model 8500 CS from Gagescope. Of course, the present invention is not 
limited to that specific embodiment of acquisition board. A typical sampling 
frequency is 500 MHz, in 8-bit format. 

The controller 16 is in the form of a personal computer 
25 including a central processing unit (CPU) 22 which is provided with an 
output device 24 in the form of a display monitor coupled to the personal 
computer 16 and input devices such as a keyboard and pointing device 
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also coupled thereto (both not shown). The controller 16 is provided with a 
memory for storing the scan signals and/or storing information elastogram 
related information as it will be explained hereinbelow in more detail. The 
controller 1 6 may take many other forms including a hand held device, an 
5 electronic circuit, a programmed chip, etc. 

The controller 16, RF signal pre-amplifier 18 and/or 
ultrasound system 11 may be part of a single vascular elastography 
device. 

10 

The controller 16 is configured and programmed so as to 
implement a method for vascular elastography as it will be described 
furthering. 



1 5 In operation of the system 1 0, the ultrasound transducer of 

the ultrasound system 1 1 is applied on the skin over the region of interest, 
and the arterial tissue is dilated by the cardiac pulsation or any other 
arterial tissue dilatation means. 



20 The elastograms, the equivalent elasticity images, are 

computed from the assessment of the vascular tissue motion as it will be 
explained hereinbelow in more detail. For that purpose, longitudinal or/and 
cross-sectional RF data are measured. For longitudinal investigations, 
because they are more convenient, axial deformation parameters may be 

25 sufficient to characterize the vessel wall. For cross-section data, the full 
strain tensor is used to compute the Von Mises parameter, because 
motion parameters might be difficult to interpret since tissue motion occurs 
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radially within the vessel wall while the ultrasound beam propagates 
axially. As a consequence of that, the elastograms are subjected to 
hardening and softening artifacts, which are to be counteracted. 

5 As it will now be described herein in more detail with 

reference to a method 100 for vascular elastography, according to a first 
illustrative embodiment of a second aspect of the present invention, the 
Von Mises (VM) coefficient is computed in order to circumvent such 
mechanical artifacts and to appropriately characterize the vessel wall. 
10 More specifically, a Lagrangian speckle model estimator (LSME) is used 
to model the vascular motion which provides the full strain tensor for 
computing the VM coefficient. 

The method 100, which is illustrated in Figures 2-3, 
1 5 comprises the following step: 

1 02 - Providing a sequence of radio-frequency (RF) images, 
including at least one pre-tissue-motion and at least one post-tissue- 
motion image, of a vessel delimited by a vascular wall; 

20 

104 - Partitioning both RF images within the vascular wall 
into corresponding data windows; 

1 06 - Approximating a trajectory between the pre- and post- 
25 tissue-motion for corresponding data windows; and 

1 08 - Using the data window trajectories to compute a strain 
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tensor in each data window. 

Each of these steps will now be described in more detail. 

In step 102, a time-sequence of one-dimensional (1D) l(x(t)), 
5 two-dimensional (2D) l(x(t), y(t)) or three-dimensional (3D) RF images 
l(x(t), y(t), z(t)) is provided, among which two images are selected for 
steps 104-108. The first image l(x(t 0 ), y(t 0 ), z(t 0 )) will be referred to as the 
pre-tissue-motion image and the second image l(x(t 0 +At), y(t 0 +At), 
z(t 0 +At)) will be referred to as the post-tissue-motion image. Images 
10 obtained through other imaging modalities than ultrasound can also be 
used. 

In step 104, both selected RF images are partitioned within 
the vascular wall into corresponding data windows Wy. Figure 4 illustrates 
15 an example of two-dimension partitioning of the region of interest (ROI) 
into W mn windows. Of course, the partitioning of the ROI can be in 1 D or 
extended in three-dimension. 

The vascular tissue and boundary conditions are generally 
20 heterogeneous. The vessel wall is thus expected to deform non-uniformly. 
As illustrated in Figure 4, to assess such tissue motion, the method 100 
includes subdividing the ROI within the vascularwall into several partitions 
Wy, for which motion can be assumed as affine. 

25 In step 106, a trajectory is approximated for each data 

windows by the zero-order and first-order terms of a Taylor-series 
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10 



expansion. For a three-dimensional tissue, assuming that the origin is set 
at (0,0,0), this can be expressed as: 
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Equation 1 defines an affine transformation, i.e. it is the 
result of a translation (vector [Tr]) and of a linear geometrical 
transformation of coordinates (matrix [LT]). Equation 1 can also be seen 
as representing trajectories that describe a tissue motion in a region of 
constant strain. Strain is usually defined in terms of the gradient of a 
displacement field; hence, as (x, y, z) represents the new position of a 
point (x 0 , yo, z 0 ), the (u, v, w) displacement vector in the Cartesian 
coordinate system is given by: 
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In this Equation, [I] is the 3D identity matrix. The ey, which are the 
components of the 3D-strain tensor e, can then be defined in terms of A 
(the deformation matrix) as: 

Bij(t) = ^MO+Ajifr)] (3) 

5 Given Equation 3, the Von Mises (VM) coefficient can be computed 
[Mase, 1970]. VM is independent of the coordinate system and is 
mathematically expressed as: 

5 = {|[( £ xx -£yy) 2 +(Syy +(^zz " ^xx f + ^ly + ^yz + *4 )]} (4) 

For each partition window within the vascular wall, £ can be 
10 related to the elastic modulus (E) through the following relationship: 



E = £ 



(5) 



where a depends on the pressure gradient that results from the blood flow 
pulsation or from the pressurization of the vessel by an external means. 



15 In step 1 08, the deformation matrix (A) is computed in each 

data window using the data window trajectories. 

More specifically, a non-linear minimization is performed for 
each Wy by computing the [LT] that allows the best match between each 
20 Wij of the pre-tissue motion image and its counterpart or corresponding 
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window in the post-tissue motion image. In other words, the method 100 
yields the deformation matrix (A) and the strain tensor (e) through 
Equations 1, 2 and 3. The map of the distribution of each component of 
the deformation matrix (A) provides a unique elastogram; the components 
5 of 8 can also be combined to provide a composite elastogram as it is the 
case for the VM coefficient (Equation 4). In Cartesian coordinates, en 
(=£xx), £22 (=£yy) and s 33 ( =e zz) are referred to as lateral, axial and transverse 
elastograms, respectively. The e 12 (=£21 =£ xy =£yx), £13 (=£31 =£xz =£zx) and 
£23 (=£32 =£yz =£zy) are the shear elastograms available with this technology. 
10 As it is believed to be well known in the art, elastograms are usually 
presented as color-code images where dark and bright regions are 
conventionally associated to hard and soft tissues. 

Step 1 08 that is referred to herein as the Lagrangian Speckle 
15 Model estimator (LSME), can be mathematically expressed as: 

MIN ||l(x(to),y(to),z(t 0 ))-I Lag (x(t 0 + At),y(t 0 + At),z(t 0 + At))f (6) 

with = [Tr;LT(:)] for a given Wy using the notation for augmented vector 
(;) and matrix vectorisation (:). Hence, vj/y is a 1 2 x 1 vector built from the 3 
x 1 Tr vector and the 9x1 vectorisation of LT. 
20 lLag(x(t 0 +At),y(t 0 +At),z(t 0 + At)) is the Lagrangian speckle image (LSI); 
it is defined as the post-tissue motion RF image 
l(x(t 0 + At),y(t 0 + At),z(t 0 + At)) thgt was numerjca || y compensated for 

tissue motion, as to achieve the best match with ^^oXy^oX^o)) 
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[Maurice and Bertrand, 1999]. The appellation "Lagrangian" refers to the 
Lagrangian description of motion. The minimum of Equation 6 is obtained 
by using the appropriate that best matches lo = l(x(t 0 ),y(to),z(t 0 )) and 
l 1 =l(x(t 0 +At),y(t 0 +At),z(t 0 +At)). v|/ 0 is the initial guess to start the iterative 
5 process. The regularized nonlinear Levenberg-Marquardt (L&M) 
minimization algorithm [Levenberg, 1963; Marquardt, 1944] is used in 
solving Equation 6. Of course, other minimization algorithms can also be 
used. 

The method 100 allows computing the full 3D-strain 
tensor (Equation 3). Whereas the divergence parameters (e xx , Syy and 8zz) 
provide information about tissue stiffness, the shear parameters (e xy , Sxz 
and 8yz) can provide useful insights on the heterogeneous nature of the 
vessel wall. 

A method 200 for vascular elastography according to a 
second illustrative embodiment of the present invention will now be 
described with reference to Figure 5. Since the method 200 is very similar 
to method 1 00, and for concision purposes, only the differences between 
the two methods will be described furthering. 

The optical flow-based method 200 is based on the 
assumption that speckle behaves as a material property. 

25 As illustrated in Figure 5, the cross-correlation analysis 

provides 3D displacement fields and a correlation map between l 0 and h. 
Tissue motion parameters (A, ty) are computed for each Wy using l Q and 



10 



15 



20 
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The material derivative of a function l(x(t),y(t),z(t)) is 



given as: 




t » 



(7) 



with I x , l y , l z and l t being the partial derivatives of l(x(t),y(t),z(t)) with 
respect to x, y, z and t, respectively. Here, l t is the time rate of change of 

l(x(t),y(t),z(t)) in the observer coordinate system, ( — j is the 



velocity vector of a "material point" located at (x,y,z), and — is the intrinsic 

1 0 rate of change of the material point. For a time-sequence of two images at 
interval dt, Equation 7 can be rewritten as: 

dl = I x dx + I y dy + I z dz + (l(x(t + dt),y(t + dt),z(t + dt)) - I(x(t), y(t),z(t)) (8) 

where (dx,dy,dz) represents the displacement vector of the "material 
point" located at (x,y,z) in the time interval dt. 

15 Furthermore, assuming speckle being a material property 

that is preserved with motion (dl(x(t), y(t),z(t)) = 0), and assuming 

l(x(t + dt), y(t 4- dt), z(t + dt)) (equivalentiy, 1^ (x(t + dt), y(t + dt), z(t + dt))) being 
a very close approximation of l(x(t),y(t),z(t)), one obtains: 
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I x dx + 1 y dy + I 2 dz = -(I Lag (x(t + dt), y(t + dt), z(t + dt)) - l(x(t), y(t), z(t))) (9) 
Now, with respect to Equations 2 and 3, Equation 9 can 

be rewritten as: 



10 



K^+A^x+A^y+A^z) (ty+A^x+Ayyy+A^z) (t z +A zx x+A zy y+A zz z)] 



5 For the purpose of clarity, the following simplifications are 

made in Equation 10: t x =x(0,0,0,t); t y =y(0,0,u,t); t z = z(0,0,0,t); and 

It = frug (*(t + dt), y(t + dt), z(t + dt)) - l(x(t), y(t), z(t))) . 

Finally, if the ROI has a size of p x q pixels, the discrete 
form of Equation 10 can be written as: 



!x, x i 



!x,yi 



I Z] X 1 



i Zl yi 



J x 2 x 2 Ix 2 Y2 Ix 2 Z2 J x 2 



J z 2 x 2 ! Z2 y2 l z 2 
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Lt 2 _ 







(11) 



I x pxq X pxq I x p ^ypxq I x p ^Vq I ^pxq "' *Vq X pxq ^Zpxq^pxq ^jxq 2 ^ ^jxq 

Solving Equation 1 1 yields A, t x , t y and t z for each pixel of the ROI. One 
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hypothesis supporting Equation 11 is that(dl(x(t),y(t),z(t))= 0), i.e. that 
the intensity of each pixel remained the same following tissue motion. In 
practice, however, because of signal decorrelation, that is rarely the case. 
Equation 1 1 then provides a solution for the minimization problem given in 
5 Equation 6. The main advantage of method 200 over the method 100 is 
relative to the processing time. Indeed, the computation time is improved 
by a factor close to 25. 

As illustrated in Figure 5, it is to be noted that this 
10 implementation of the LSME uses cross-correlation analysis to compute 
motion compensation as to provide I Lag (x(t + dt),y(t + dt),z(t + dt)). 

Interestingly, a total of 16 parameters can be assessed in 3D 
(9 parameters in 2D). Nevertheless, the strain parameters are so far the 
most convenient for the purpose of characterizing soft tissue mechanical 
15 properties. 

Non-Invasive Vascular Elastoqraphv (NIVE and MicroNive) 

Tissue motion analysis for cross-section data 

20 

Superficial arteries such as the carotid and femoral are easily 
accessible and can be imaged longitudinally. This can be seen as the 
most convenient application of the method 100 in non-invasive vascular 
elastography (NIVE), since tissue motion can be expected to run close to 
25 parallel to the ultrasound beam. In this context, the axial components of 
the deformation matrix may be sufficient to characterize the vessel wall. 
Whereas tissue motion analysis for longitudinal data will be presented with 
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reference to a further illustrative embodiment of the present method, we 
here emphasize on cross-sectional data. 

For a continuum, motion can be described in a Lagrangian 
5 coordinate system or in an Eulerian coordinate system [Le Mehaute, 1 976; 
Hinze, 1975]. In the literature, the Eulerian coordinate system is 
sometimes referred to as the observer's coordinate system, whereas the 
Lagrangian coordinate system is known as the material coordinate 
system. The material coordinates allow expressing each portion of the 
1 0 continuum as a function of time and position. Whereas the mathematical 
formulation (Equations 1-4) is in 3D, the illustration is here presented in 2D 
for the purpose of simplification. The difference between these two 
coordinate systems is illustrated in Figure 4, where the (x,y) constitutes the 
observer's coordinates and the (r,cp) defines the material coordinates. 

15 

With most imaging systems, such as ultrasonography, the 
observer's and the material coordinate systems are generally the same; 
hence, most tissue motion estimators use, by definition, the observer's 
coordinate system. However, the material coordinates can be presented 
20 as a suitable way to describe speckle dynamics [Maurice and Bertrand, 
(1999)]. 

As illustrated in Figure 4, the observer's coordinate system is 
the Cartesian (x,y)-plane. This system is different from the motion 
25 coordinate system that is in the radial (r,cp)-plane. In such a situation, the 
parameters of an estimator are expected to be very difficult to interpret. 
One of the challenges of non-invasive vascular elastography, regarding 
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5 



cross-section data, concerns the interpretation of the estimated motion 
parameters to characterize the vascular tissue. 

Simulated and Experimental in vitro Results of NIVE and MicroNive 

As it will now be described in more detail with reference to a 
pathology free simulation, the VM coefficient can be used as a tissue 
characterization parameter to better interpret the displayed images. 

Motion analysis for a homogenous tissue 



A pathology-free application simulation will now be 
10 considered, that is the case of a circular, axis-symmetric and 
homogeneous vessel section. To take into account the constraints 
induced by the environmental tissues and organs, it is hypothesized that 
the vessel section is embedded in an infinite medium of higher Young's 
modulus. A pressurized thick-wall cylindrical blood vessel of inner and 
15 outer radii Rj and R Q respectively, embedded in an elastic coaxial 
cylindrical medium of radius R e , is considered. It is assumed that the plane 
strain condition for the vessel wall applies and also that the two media are 
incompressible and isotropic. Referring to Equation 2, the displacement 
gradient components (Equation 2) are given by: 



A(x,y)= 



20 



K -2K Xy 



(x 2 + y 2 ) 2 "(f + y 2 ) 2 

"(x'H-y*)' •(.'*,')' j 



(12) 
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with: 



K„ =-R 



R 




(13) 



where the subscripts (1 ) and (2) describe respectively the vessel wall and 
the external medium; Pb defines the blood pressure and E is the Young's 
modulus. 



Equations 12 and 13 were implemented to simulate the 
dynamics of a homogeneous vessel section subjected to an intraluminal 
pressure. Close to 6 % intraluminal dilation was induced regarding the 
10 constitutive model presented here. It is to be noted that 6 % intraluminal 
dilation is equivalent to 3 % compression of the intraluminal wall. The 
physical vessel dimensions were 7-mm outer diameter and 4-mm inner 
diameter as to approximate the physiological case of a femoral artery. 

Figures 6A and 6B present respectively the lateral and axial 
displacement fields; they include gray-scale "colorbar" expressing the 
displacement in jam (10' 6 m). Maximum motion occurred at the lumen 
interface. Figures 6C to 6F present the Ay components of Equation 12, 
which are respectively the lateral strain, the lateral shear, the axial shear 
and the axial strain; they include gray-scale "colorbar" expressing the 
strain in percentage. 



15 



20 



Ayy is expected to be less or equal to zero (< 0) since, in 
conventional elastography, an external force is applied and induces tissue 
25 compression. Traditionally, smaller strain amplitude values are associated 
with harder regions and are printed dark; equivalently, higher strain 
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amplitude values are associated with softer regions and are printed bright. 
However according to the method 1 00, as it can be observed in Figure 6F, 
dilation can also be detected (Ayy > 0) in the elastogram. In an 
elastographic sense, the dilation regions can be misinterpreted as soft 
5 tissue. Indeed, in Figure 6F, two harder zones (Ayy < 0) likely seem to be 
identified at 12 and 6 o'clock. Because, for the conditions simulated, the 
vessel wall is homogeneous; such a phenomenon is referred to as 
hardening artifact [Ophir et a/., (1999)]. Inversely, two softer zones (Ayy > 
0) at 3 and 9 o'clock seem also to be identified; they are referred to as 
10 softening artifacts. Such motion artifacts stem from the fact that motion 
occurs radially and is observed in Cartesian coordinates. 

The motion parameters in their natural polar coordinate 
system or material coordinate system will now be considered. In Figure 

15 7A, the radial displacement field is computed from the lateral and axial 
displacement fields (Figures 6A and 6B respectively). The radial 
displacement field is also presented in a polar (r,cp) coordinate system 
(Figure 7B). The gradient of the latter displacement field thus provides the 
radial strain (Figure 7C). Figure 7D shows a plot of the radial strain at cp = 

20 7i. One can observe the monotonic profile of this plot, being maximal at the 
lumen and minimal at the outer side of the vessel. Such a phenomenon is 
a consequence of the geometry and is known in the literature as the strain 
decay [Ryan and Foster, (1997)]. Finally in Figure 7E, the radial strain is 
reported back in the (x,y) coordinate system. Regardless of the strain 

25 decay phenomenon, one can observe that no specific hard or soft region 
is identified. Figure 7E thus illustrates a strain profile that adequately 
represents a homogenous vessel wall behavior. 
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Hopefully, elastograms such as the one shown in Figure 7E 
allows to appropriately characterizing the vessel wall. However, in NIVE or 
MicroNive, motion is studied in the transducer coordinate system; that is 
5 the (x,y)-Cartesian coordinates. Accordingly, elastograms are expected to 
be as artifactual as the one in Figure 6F. The VM coefficient (Equation 4) 
is then used to characterize the vessel wall [Mase, 1970]. 

A comparison between the radial strain and the Von Mises 
10 parameter (^) is shown in Figures 8A-8B for a homogeneous vessel wall. 
Qualitatively, both parameters are equivalent. Figure 8C shows the plots 
for the radial strain ( — ) and % ( — ) at x = 0. Those graphs show that the 
VM coefficient likely improves the contrast between higher and lower 
strains while the profile remains the same. Moreover, regardless of the 
15 strain decay, £, (as well as the radial strain) is interestingly free of 
hardening or softening artifacts and is thus suitable to non-invasively 
characterize the vessel wall. To corroborate such an assumption, a more 
complex geometry, that is a heterogeneous vessel wall, has been 
investigated in vitro and the results are presented in the following section 

20 

MicroNIVE experiments have been conducted on vessel- 
mimicking phantoms with a lumen diameter of 1.5 mm and a wall 
thickness of 2 mm. The phantoms were made of polyvinyl alcohol cryogel 
(PVA-C). 

25 

An experimental set-up 26, used to produce mechanical 
deformation of polyvinyl alcohol cryogel (PVA-C) vessel-mimicking 
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phantoms, and including the system 10 to collect RF ultrasound data that 
can be used in computing vascular elastograms according to the method 
from 100 will now be described with reference to Figure 9. 

5 A mixture of water-glycerol was circulated in a flow phantom 

30. The height difference between the top and bottom reservoirs 28 and 
36 allowed adjustment of the gravity-driven flow rate and static pressure 
within the lumen of the phantom 30. A peristaltic pump 38 was used to 
circulate the fluid from the bottom to the top reservoirs 36 and 28. The flow 

10 rate was measured with an electromagnetic flowmeter 32, which was a 
Cliniflow II, model FM 701 D from Carolina Medical, and the pressure was 
monitored by a MDE Escort instrument 34, which was a model E102 from 
Medical Data Electronics. As illustrated in Figure 9, the flow phantom 30 
was not directly connected to the tubing of the top reservoir 28 to facilitate 

15 the small incremental pressure step adjustments necessary to obtain 
correlated deformation of the RF signals within the PVA-C vessel wall. 

As can be seen in Figure 10, the polyvinyl alcohol cryogel 
PVA-C vessel 39 of the flow phantom 30 was positioned between two 
20 watertight connectors 40, in a Plexiglas box 42 filled with degassed water 
44 at room temperature. Rubber o-rings were used to tight the PVA-C 
vessel 39 onto Plexiglas tubes 46 at both extremities. 

Based on previous works by Chu and Rutt (1997), the tissue- 
25 mimicking vessel 39 was made of PVA-C. This biogel solidifies and 
acquires its mechanical rigidity by increasing the number of freeze/thaw 
cycles. Indeed, the number of freeze/thaw cycles modifies the structure of 
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the material by increasing the reticulation of fibers. It has been shown that 
the elastic and acoustic properties of PVA-C are in the range of values 
found for soft biological tissues [Chu and Rutt, 1997]. More specifically, it 
has been demonstrated that the stress-strain relationship can be very 
5 close to that of a pig aorta. 

The vessel-mimicking phantoms 30 approximately had a 
1 .5-mm lumen diameter, 2-mm wall thickness, and 52-mm length. A 1 .5% 
by weight of Sigmacell (type 20, #S-3504, from Sigma-Aldrich) was added 
10 to the PVA-C to provide acoustical scatterers. 

Results for one double-layer vessel will now be presented. 
Each layer had a thickness close to 1 mm, and the inner layer was made 
softer than the outer one. The numbers of f reeze-thaw cycles were set at 
15 2 and 4 for the inner and the outer portions of the wall, respectively. Each 
freeze-thaw cycle took 24 hours and the temperature was incrementally 
. varied from -20 C to 20 C, by using a specifically designed electronic 
controller (Watlow, model 981) and a freezer equipped with heated 
elements such as Supra Scientifique's model YF-204017. 

20 

Figures 11A-11C show a schematic representation of the 
moulds that were used to construct the double-layer vessel-mimicking 
phantoms 39, the simulated vessel having a 1 .5 mm lumen diameter, a 2 
mm wall thickness (roughly 1 mm for each layer), and a 52 mm length. At 
25 a first instance, PVA-C was poured between the first and second 
templates; that underwent (n Q - r\\) freeze/thaw cycles to provide the 
external layer. For the purpose of clarity, n\ and n Q as the numbers of 
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cycles for the inner and the outer layers, respectively. Secondly, fresh 
PVA-C was poured between the second and third templates, while 
maintaining the first template in place; that underwent nj freeze/thaw 
cycles to provide a complete double-layer vessel-mimicking phantom. 

5 

Returning to Figure 9, the ultrasound biomicroscope 12 
(Visualsonics, model VS-40) provides an RF output from which the 
received RF data were transferred to a pre-amplifier 18 (Panametrics, 
model 5900 PR). After amplification, the signals were digitized with an 
10 acquisition board 14 (Gagescope, model 8500 CS) installed in a personal 
computer 12. The sampling frequency was 500 MHz, in 8-bit format. 

The double layer vessel-mimicking phantoms 30 measured 
5.5 mm in outer diameter, whereas the RF images extended to 8 mm x 
15 8 mm. Measurement-windows (partitions or ROI) were of 272 jam x 
312 jam (200 samples x 20 RF lines), with 85 % axial and lateral overlaps. 
The estimated motion parameters were post-processed using a 5*5 kernel 
Gaussian-filter. The pressure pre-load was 10 mmHg, and the pressure 
gradient was 5 mmHg between subsequent images. 

20 

Figure 12A shows a B-mode image (at 10 mmHg) of the 
phantom 39. To ensure the reproducibility of the method, VM elastograms 
(Equation 4) were computed comparing im^ and im 2 j (j = 1 , 2, 3, 4). Figure 
1 2B presents an average of 4 such elastograms that shows the visibility of 
25 both layers. Figure 12C shows a graph of an average of 5 axial lines 
chosen in the middle of the VM elastogram (around x = 0). Both layers can 
quantitatively be differentiated, specifically at the upper side of the vessel 
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(left side of panel c) where maximum strains around 0.8% and 0.4% are 
observed for the inner and the outer layers, respectively. 

As shown by Equation 4, the composite VM elastogram is 
5 obtained by computing the four components of the deformation matrix. 
Whereas the method 100 adapted to NIVE or MicroNive provides very 
accurate axial deformation estimates, lateral motion assessment remains 
a non-trivial problem, due to the relatively limited lateral resolution of RF 
ultrasound images. This difficulty was partially circumvented. Indeed, 
10 under the assumptions of weak compressibility of biological soft tissues 
and isotropy inside each measurement-window, the lateral strain (e^ = Axx) 
was deduced from the axial strain (syy = Ayy), that is Ayy = Axx/2. 

From Figure 12B, the strain in the inner layer close to the 
15 vessel lumen was on average 1.11 ± 0.05 %. Since the intraluminal 
pressure gradient was 5 mmHg, the elastic modulus E was estimated at 
60 ± 3 kPa for this material (made with two freeze-thaw cycles). It is to be 
noted that the elastic modulus E, for the inner layer, was estimated from 
Equation 5. Indeed, as a first approximation, a for this layer is given by the 
20 static pressure gradient inside the vessel measured for the conditions 
corresponding to the pre-motion and post-motion RF images. E has been 
estimated at around 49 ± 6 kPa for a 1 freeze-thaw cycle PVA-C. In both 
cases, the pressure pre-load was 10 mmHg. The elastic modulus E was 
higher for the 2 freeze-thaw cycles material as it could be expected, since 
25 PVA-C made of 1 freeze-thaw cycle is softer than PVA-C made of 2 
freeze-thaw cycles. 
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Even though the experimental set-up from Figure 9 has been 
described including a high-resolution transducer, with a 32 MHz central 
frequency, to acquire the RF data for the in vitro experimentations of 
MicroNive, lower frequencies can also be used, for example 10 MHz or 
5 lower for in vivo NIVE applications to counteract ultrasound signal 
attenuation. Therefore a transducer is carefully selected for the system 10 
since it is considered that, for ultrasound signals, lower is the frequency, 
deeper the beam propagates in the tissue, but, higher is the frequency, 
better is the resolution. 

10 

NIVE applications of the method 100 and of the system 10 
include characterizing abdominal or peripheral aneurysms and superficial 
arteries such as the femoral and the carotid. 

15 In vivo experimentation of NIVE on the carotid artery of a normal 
human subject 

RF data were acquired from the carotid artery of a young 
healthy human subject. The Ultrasonix™ ES500RP ultrasound system, 

20 along with a 7.5 MHz transducer, was used to image longitudinal sections 
of the artery. Figure 1 3A shows a B-mode image of the artery. Figure 1 3B 
presents a manual segmentation of the same artery. Figures 13C-13D 
present axial strain elastograms computed with the method 100; the gray- 
scaled "colorbar" expresses the deformation in percentage. Since these 

25 elastograms were computed from data acquired during diastole, the axial 
strain values are expected to be positive. The regions of interest 
highlighted in Figures 13C and 13D correspond to sections of the carotids 
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where motion occurred close to parallel to the ultrasound beam. 
Interestingly, in both elastograms, the upper vessel walls are observed to 
deform less than the lower walls; that is because the force exerted by the 
transducer over the skin can be seen as a boundary condition limiting the 
5 motion of the upper vascular tissues. It is to be noted that the Von Mises 
coefficient has not been used to display the strain patterns obtained from 
the method 100, because longitudinal sections of the carotid vessels were 
acquired instead of transverse planes. 

10 In vivo experimentation of MicroNIVE on the carotid artery of 
normotensive and hypertensive rats 

The method 100 can be used to characterize mechanical 
properties of small vessels (MicroNIVE) in humans or animals. More 
1 5 specifically, the method 1 00 can be used in the context of the phenotyping 
in hypertension (HT) with genetically-engineered rat models. 

High-frequency ultrasound RF data were acquired on 6 male 
rats: 3 normotensive Norway Brun rats (labeled as NT1, NT2 and NT3) 

20 and 3 spontaneously hypertensive SHR rats (HT1, HT2 and HT3), 
respectively. All animals were 15-weeks old and they were anesthetized 
by inhalation of 1.5 % isofluorane during RF data acquisition. The body 
temperature of each animal was monitored with a rectal probe and 
maintained at 37 ± 1 °C by using a heating surface. The hairs over the 

25 neck were shaved and further removed with a depilatory cream. 
Sequences of RF data were acquired over longitudinal sections of the 
carotid artery using the Vevo 660™ from Visualsonics, an improved 
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version of the VS-40 ultrasound system (that was used to acquire the in 
vitro phantom data reported above in Figure 12). Indeed, this ultrasound 
system is equipped with an encapsulated oscillating element transducer 
and provides a frame rate of 30 images/s. For the purpose of the in vivo 
5 investigation, a 40 MHz central frequency transducer was used. The mean 
systolic pressures for NT and HT rats (measured with a tail-cuff monitoring 
system) were 87 ± 12 mmHg and 158 ± 16 mmHg, whereas the heart 
beats were 323 ±15 beats / min. and 365 ± 6 beats / min., respectively. 

1 0 Elastograms were computed using the method 1 00 adapted 

for MicroNIVE as it will now be described. All successive acquired RF 
images that were digitized over several cardiac cycles were used. No 
averaging was used to display the axial elastograms of Figures 14C-14H. 
Manual segmentation has been done to display only the strain patterns 

15 within the vascular wall. It is to note that the Von Mises coefficient has not 
been used to display the strain patterns obtained from the method 100, 
because longitudinal sections of the carotid vessels were acquired instead 
of transverse planes. 

20 Figures 14A-14B show two B-mode images obtained for a 

normotensive rat (NT1) and a hypertensive one (HT2), respectively. In 
both cases, the internal diameter of the carotid was around 1.1 mm, 
whereas the wall thickness was close to 160 jim for NT1 and 120 jim for 
HT2, respectively. 

25 

Figures 14C-14H show axial strain elastograms computed 
using the method 100; the gray-scale "colorbar" providing the strain in 
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percent. The negative strains are indicative of vessel dilation (diastolic 
phase). Since it has been difficult, for most rats, to have longitudinal 
sections of 6 mm, only portions of the carotids are displayed on the 
elastograms. The carotids of the three normotensive rats (NT1 , NT2 and 
5 NT3) appear on average twice softer (strain values up to 7 %) than those 
of the hypertensive ones (HT1 , HT2 and HT3), where a maximum of 3.3 % 
strain was estimated. To provide a more rigorous interpretation of these 
data, one should of course consider the pressure gradient between the 
pre-compression (pre-tissue-motion) and post-compression (post-tissue- 
10 moition) RF images used to compute the elastograms. This pressure 
gradient was higher for HT than NT rats. However the exact values are not 
known since the RF images were acquired without any ECG or pressure 
synchronization.. 



A method and system for MicroNIVE according to the 
present invention can be used in ex-vivo experiments or in vivo testing on 
animals or humans. For example, using recombinant inbred strain (RIS) 
rats, the method can be used to examine the modulation of drug-induced 
cardiovascular remodeling as a function of HT and aging. Examples of 
protocols for ex-vivo and in vivo experiments are described in the 
following. 

Animals are treated with placebo, losartan, which is an 
antihypertensive drug and an antagonist of angiotensin II (ANG II) type 1 
25 (AT1) receptors (30mg/day), and nifedipine, an antihypertensive drug, 
which is a calcium channel blocker (30 mg/day) for two weeks starting at 
12 weeks of age. Regarding the ex-vivo experiments, the animals are 



15 



20 



WO 2005/074804 



PCT/CA2005/000162 



41 



killed and segments of arteries (carotid, for example) are excised. 
Segments (« 2-cm in length) will be mounted on similar apparatus than for 
the vessel-mimicking phantom experimentation described above in 
Figures 9 and 10. Examples of preparation and of protocol of study for 
5 such arteries are well documented [Li et al., 1998 and 2003; Intengan et 
al., (1998a and 1998b)]. The vessel is adjusted to its length before 
excision such as the vessel walls become parallel. The vessel is 
equilibrated under a constant intraluminal pressure of 45 mmHg with 
physiological salt solution [Intengan et al., (1998a and 1998b)]. A 

10 servocontrolled pump stepwise (5 mmHg each step) increases the 
intraluminal pressure, and time-sequence RF data are acquired at 
different frequencies (25 or 40 MHz, depending on the artery) with an 
ultrasound biomicroscanning system, such as the Vevo 660TM from 
Visualsonics. The elastograms are computed using a method for vascular 

1 5 elastography according to the present invention, such as the method 1 00. 

Other In vivo experiments can also be performed using RIS 
rats. These animals are treated with placebo, losartan (30mg/day), and 
nifedipine (30 mg/day) for two weeks starting at 12 weeks of age for the 

20 purpose of examining the modulation of drug-induced cardiovascular 
remodeling as a function of HT and aging. To acquire RF data, the rats 
are anesthetized by inhalation with 1.5 % isofluorane. Physiological 
parameters (temperature, pressure and ECG) are monitored. The 
temperature is maintained close to 37°C using a hot plaque. The region of 

25 interest is shaved using a conventional electric shaver; the remaining hair 
is removed with NairTM or another lotion hair remover. Even though there 
is possible impact of anesthesia on the cardiac function of rats (reduction 



WO 2005/074804 



PCT/CA2005/000162 



42 



of the heart beat, cardiac output, etc.), and possibly on the arteries, that 
effects is the same for all animals and therefore does not interfere with the 
interpretation of the results since all the rat strains are anesthetized. 

5 In the case where carotid arteries are investigated, only the 

axial strain are required to study longitudinal images, since the vessel wall 
motion can be seen as running parallel with the ultrasound beam. The RF 
data are processed using the method 100 to provide step-wise 
elastograms (strain images). From the strain estimates, another 
10 mechanical parameter (namely stress/strain ratio) is calculated. 

The MicroNIVE method according to the present invention 
allows providing significant new insights regarding the pathophysiology of 
HT and aims at leading to new discoveries in the field of pharmacology for 
15 example, even though it is not limited to this particular application. 

Simulated and Experimental in vitro Results on Endovascular 
Elastography (EVE) 

20 A method for endovascular elastography (EVE) according to 

a third illustrative embodiment of the second aspect of the present 
invention will now be described. Since the method according to this third 
embodiment is similar to the method illustrated in Figure 2, and for 
concision purposes, only the differences between these two methods will 

25 be described herein. 

The first step of the method is to acquire intravascular RF 
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images using a catheter. Following the example of IVUS, and as 
schematically illustrated in Figure 15, a transducer is placed at the tip of 
the catheter and cross-sectional imaging of a vessel is generated by 
sequentially sweeping the ultrasound beam over a 360° angle. It is to be 
5 noted that, in the ideal situation illustrated in Figure 15, the ultrasound 
beam runs parallel with the vascular tissue motion, i.e. in the (r,<p) 
coordinate system. 

Mechanical parameters (radial strain, in this case) are then 
1 0 estimated from analyzing the kinematics of the vascular tissue during the 
cardiac cycle or in response to an angioplasty-balloon push or to any other 
force exerted axially onto the inner vascular wall. 

The Lagrangian Speckle Model Estimator (LSME) is then 
1 5 formulated for investigations in EVE, i.e. using a polar coordinate system. 
Indeed, while the full 2D-deformation matrix A. can be assessed, only the 
radial strain component A rr (= e rr ) is displayed. This is motivated by the fact 
that tissue motion, in EVE, is expected to run close to parallel with the 
ultrasound beam. Again, the LSME allows computing the full deformation 
20 matrix A of Equation 3. Since A is directly assessed, no derivative of the 
displacement fields is required, as it was also the case for the first and 
second embodiments of the present invention. Although the method is 
general and applies to either 1 D, 2D or 3D RF data, the description given 
hereafter refers to a 2D model for simplicity. 

25 

In 2D, using polar coordinates, the LSME can be formulated 

as: 



WO 2005/074804 



PCT/CA2005/000162 



44 



MIN 



l(r,cp,0)-[l(r 5 <p,t)] LT f = MIN |l(r,(p ; 0)-I Lag (r,(p,t) 

P II T T 1 



LT p 

= MIN| M ( r ^ t 



.(14) 



The minimum is obtained by using the appropriate [LT P ]. 
[LTp] is a linear transformation matrix which maps the Cartesian 
trajectories in a polar coordinate system. However, for a small ROI (Ar, Aq>) 
5 that is far from the vessel lumen center, motion equivalently can be 
investigated using either a polar or a Cartesian coordinate system. In other 
words, the following approximation can be done to compute the 
elastogram: 

^ = LT-I = LT p -l (15) 

10 where I is the 2D-identity matrix. In other words, the solution to Equation 
13 can be obtained from solving Equation 6. 



Biomechanical simulations of vessel wall kinematics 



A computational structural analysis has been performed on 
15 one simulated idealized coronary plaque (see Figure 16) and on a model 
identified on Figure 17B created from measurements made of a typical 
composite plaque identified from an in vivo IVUS image of a patient with 
coronary artery disease (see Figure 17A). The former allowed validating 
the potential of the EVE method according to the present invention to 
20 differentiate between hard and soft vascular tissues and the latter allowed 
characterizing the heterogeneous nature of atherosclerotic plaques, which 
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is linked to the risk of rupture and thrombosis. 

For the two models, the materials were considered as quasi- 
incompressible (Poisson ratios v = 0.49) and isotropic with linear elastic 
5 properties. The Young's modulus for the healthy vascular tissue (or 
adventitia & media) was 80 kPa [Williamson et a/., (2003)], while the 
dense fibrosis (much stiffer) was set at 240 kPa, and the cellular fibrosis 
(softer than the dense fibrosis) was chosen at 24 kPa [Ohayon et al., 
(2001); Treyve et al., (2003)]. Whereas the surrounding tissue was not 
10 investigated, the bulk boundary conditions, as it may eventually be 
provided by surrounding organs, were simulated by imbedding the vessel 
in a stiffer environment of 1000 kPa Young's modulus. 

Finite element (FE) computations were performed by 
15 considering static simulations of coronary plaques under loading blood 
pressure. The simulations were performed on the geometrical models 
previously described (see Figures 16 and 17B). Nodal displacements were 
set to zero on the external boundaries of the surrounding tissue. The 
various regions of the plaque components were then automatically 
20 meshed with triangular (6 nodes) and quadrangular (8 nodes) elements. 
The FE models were solved under the assumption of plane and of finite 
strains. The assumption of plane strain has been made because axial 
stenosis dimensions were of at least the same order of magnitude as the 
radial dimensions of the vessel. Moreover, the assumption of finite 
25 deformation was required as the strain maps showed values up to 30% for 
physiological pressures [Loree etal., (1992); Cheng et a/., (1993); Lee et 
al., (1993); Ohayon et a/., (2001); Williamson etal., (2003)]. However, the 
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kinetics reported were achieved with small pressure gradients (around 
15 mmHg) such that the radial strain remained below 10 %. The Newton- 
Raphson iterative method with a residual nodal tolerance of 4x1 0' 4 N was 
used to solve the FE models. The calculations were performed with a 
5 number of elements close to 7200. This computational structural FE 
analysis was used to perform the kinematics of the vascular tissue. 

Adynamic image-formation model was implemented using 
the Matlab™ software. It is assumed that the image formation can be 
10 modeled as a linear space-invariant operation on a scattering function, 
and that the motion occurs in plane strain conditions (such as no 
transverse deformation is involved). It considers a scattering function 
Z(x(t 0 ),y(t 0 )) (at t 0 =0) that simulates the acoustical characteristics of a 

transverse vascular section in Cartesian coordinates. Knowing the axial 
15 and lateral displacement fields, they are applied upon Z(x(t 0 ),y(t 0 )) to 
perform tissue motion, thus providing Z(x(t),y(t)). The last step consists of 
convolving Z(x(t),y(t)) with the PSF (point-spread-function) to provide a 
dynamic sequence of RF images l(x(t),y(t)) or equivalents l(x,y ? t). The 
PSF is the equivalent image of a single cellular ultrasound scatterer. In 
20 other words, the PSF expresses the intrinsic characteristics of the 
ultrasound imaging system. It can be determined experimentally by using 
a phantom (a box containing a tissue-mimicking gel) containing a point 
target. The dynamic image-formation is of interest to simulate the RF 
data. 

25 

The idealized vessel illustrated in Figure 16 measured about 
3.8 mm in outer diameter, whereas the RF images extended to 4 mm x 
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4 mm. The real case vessel illustrated in Figure 17B measured about 

7 mm in outer diameter, whereas the RF images extended to 8 mm x 

8 mm. For the purpose of simulations, the intraluminal pressure gradients 
were set at 15.79 mmHg and 11.73 mmHg for the idealized and the 

5 realistic vessels, respectively. Considering the above, the dilation at the 
inner wall was around 7 % in both cases. The PSF characterized a 20 
MHz central frequency IVUS transducer. The LSME was implemented to 
assess tissue motion. Measurement-windows of 0.38 mm x 0.40 mm and 
0.77 mm x 0.80 mm, with 90 % axial and lateral overlaps, were used for 
10 the idealized and the realistic cases, respectively. 

Investigation of the "ideal" plaque pathology 

Figure 18A presents the theoretical radial strain elastogram, 
computed for the "ideal" pathology case, using Ansys FE and Matlab 

15 softwares. The plaque can slightly be differentiated from the normal 
vascular tissue, whereas a region of higher strain values is observed at the 
right portion of the inner vessel wall. This "mechanical artifact" is a direct 
consequence of the well known strain decay phenomenon [Shapo et a/., 
(1996a)]. For a more quantitative illustration, are presented in Figure 18B 

20 plots from the theoretical elastogram for two orthogonal orientations along 
x and y. Indeed, the vertical plot ( — ) shows low contrast between the 
plaque and the normal vascular tissue, whereas the horizontal plot ( — ) 
clearly points out the presence of strain decay. 

25 Figure 18C presents the radial strain elastogram as 

computed using the EVE method from the present invention, using 
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simulated RF images. As for the theoretical elastogram in Figure 18A, the 
plaque is slightly distinguishable from the normal vascular tissue. The 
graphs of Figure 17D confirm such an observation. 

5 As illustrated in Figures 18A and 18C, the present 

invention allows both characterizing the strain in the vessel quantitatively 
in addition to qualitatively. Indeed, the gray-scale "colorbars" at the right of 
each Figure express the strain in percent. 

For the purpose of compensating for strain decay, the 
1 0 radial strain elastogram resulting from the method according to the present 
invention was post-processed. Indeed, e rr was modulated with a function 
proportional to the square of the vessel radius. The strain-decay- 
compensated elastogram issued from the EVE method according to the 
present invention is represented in Figure 19A and shows substantial 
1 5 contrast improvement. For instance, the axial plot of Figure 1 9B shows an 
effective contrast ratio close to 3 between the plaque and the normal 
vascular tissue, as it can be expected. Equivalently, Figure 19C also 
shows some valuable contrast ratio improvement compared to Figure 18D. 

20 Investigation of "realistic" vessel wall pathology 

Figure 20A illustrates the theoretical radial strain elastogram, 
computed for the "realistic" pathology case. Interestingly, complex strain 
patterns are observed; nevertheless, different regions can be identified. 
25 For instance, since the ratio of Young's moduli between the dense and the 
cellular fibroses was set to 10, both of those materials can be 
distinguished. Less contrast is seen between the cellular fibrosis and the 
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healthy vascular tissue because their Young's modulus contrast was set 
to 3. As illustrated with vertical and horizontal 1 D plots from the 
elastogram (Figures 20B and 20C, respectively), strong strain decay is 
observed specifically at the inner portion of the vessel wall. 

5 

Figure 21 A illustrates the radial strain elastogram as 
computed using the method for endovascular elastography according to 
the third illustrative embodiment of the second aspect of the present 
invention, using simulated RF images. Comparing to the theoretical 

1 0 elastogram in Figure 20A, very complex strain patterns are also observed. 
Moreover, the dense and the cellular fibrosis tissues can be identified. 
However, while less prominent than in the "ideal" case study, strain decay 
remains a significant factor to compensate for to improve image 
interpretation. This is illustrated in Figures 21 B and 21 C, where vertical 

15 and horizontal 1 D graphs from the elastogram are presented. 

Figure 22A illustrates the strain-decay-compensated LSME 
elastogram, showing substantial contrast improvement. Both the vertical 
graph (Figure 22B) and the horizontal one (Figure 22C) show more 
20 effective contrast ratio between dense and cellular fibroses, and between 
cellular fibrosis and the normal vascular tissue. Moreover, it is interesting 
to notice the presence of moderate strain values (around 0.6 to 0.8 %) at 
the extremities of the plots; this characterizes regions of healthy vascular 
tissue, namely the media and adventitia. 

25 

The method for endovascular elastography according to the 
present invention has also been validated in vitro using a fresh excised 
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human carotid artery. The experimental set-up 50 used in the validation is 
illustrated in Figure 23. The set-up 50 includes a system 52 for 
endovascularelastography according to a second embodiment of the first 
aspect of the present invention. 

5 

The system 52 comprises an ultrasound scanner 54 in the 
form of a CVIS (ClearView, Cardiovascular Imaging System Inc.) 
ultrasound scanner, working with a 30 MHz mechanical rotating 
single-element transducer (not shown), a digital oscilloscope 56, more 
1 0 specifically the model 9374L from LECROY, and a pressuring system 58. 

The extremities 60-62 of an artery 64 are fixed to two rigid 
sheaths by watertight connectors 66, separated according to the original 
longitudinal dimension of the vessel 64 before excision. The intravascular 

15 catheter 68, part of the system 52 was introduced through the proximal 
sheath into the lumen of the artery 64, and then through the distal sheath. 
The distal sheath was closed with a clamp 70 to insure watertightness of 
the system 58. Injecting fluid inside the system 58 resulted in an increase 
of the pressure inside the arterial lumen since the sheath is rigid and the 

20 system is watertight. A syringe 72 is then connected to the proximal 
sheath and the inner pressure is increased or decreased by manually 
varying the fluid volume (precision: AV= 0.01 ml) inside the lumen. 
Whereas the quantitative pressure values were not monitored by an 
independent means, the fluid volume inside the lumen was maintained 

25 constant during each acquisition of RF ultrasound data. 



The ultrasound probe 74 was fixed approximately at the 
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center of the arterial lumen by two guiding elements. This protocol was 
used to limit probe motion and accordingly to reduce geometrical artifacts 
[Delachartre et al. (1999)]. 

5 A sequence of radio-frequency (RF) images was collected 

while incrementally adjusting the intraluminal static pressure steps. At 
each static pressure step (volume step), a scan of 256 angles was 
performed. A set of 11 RF images was so acquired for consecutive 
increasing physiologic fluid pressure levels. Sampling of the data was 

10 phase-synchronized, with the top image synchronizer and the RF signal 
synchronization (external outputs of the CVIS ultrasound scanner). The 
top image synchronizer allows the user to select an angular position from 
which the acquisition started; it thus permitted the acquisition of sets of 
images angularly aligned. The RF signal synchronization was done at the 

1 5 pulse repetition frequency of the bursts transmitted to the single-element 
transducer. RF data were digitized at a 500 MHz sampling frequency in 
8 bits format, stored on a PCMCIA hard disc in the LeCroy oscilloscope 
and processed off line. 

20 As shown by histology (Figures 24A and 24B), the artery was 

characterized by a thin atherosclerotic plaque (located at about 3 o'clock), 
that was only restricted to a confined angular sector. The coloration with 
saffron haematoxylin-eosin revealed that the plaque contained cholesterol 
crystals and inflammatory cells. It is to be noted that the IVUS image on 

25 Figure 24C does not clearly allow differentiating the plaque from the 
healthy vascular tissue and therefore appears insufficient to characterize 
vascular tissue. 
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Figures 25A-25J show 10 radial elastograms that were N 
computed, using the set of 11 RF images acquired for consecutive 
increasing physiologic fluid pressure levels using the method for 
5 endovascular elastography according to the present invention. 



The elastogram obtained for the lowest intraluminal pressure 
(i.e. from the 1 st and 2 nd RF images, in this case) is displayed in Figure 
25A, whereas Figure 25J shows the elastogram for the highest pressure 

10 difference (i.e. the elastogram computed with the 1 st and 1 1 th RF images). 
Indeed, maximum strain values close to 0.6 % are observed in Figure 25A, 
whereas the maximum is close to 3 % in Figure 25J. To summarize, 
elastograms in Figures 25A and 25J are the least representative, and 
those from Figure 25C to Figure 25E present very good plaque 

15 detectability, accuracy in plaque dimensions, and significant contrast 
between plaque and surrounding tissue. This demonstrates that a range of 
intraluminal pressures for which tissue motion estimation appears optimal 
exists. 

20 It is to be noted that, in the elastograms from Figures 25A- 

25J, lateral and axial values are dimensions in centimeters, while the gray- 
scaled "colorbars" give the strain in percent. 



The above results showed the potential of the method for 
25 endovascular elastography according to the present invention to 
characterize and to distinguish an atherosclerotic plaque from the normal 
vascular tissue. Namely, the geometry as well as some mechanical 
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characteristics of the detected plaque are in good agreement with 
histology. The results also suggested that there might exist a range of 
intraluminal pressures for which plaque detectability is optimal. The plaque 
at 3 o'clock displayed low strain values indicative of a stiff tissue. 

5 

The EVE method according to the present invention further 
allows providing quantitative parameters to support clinicians in diagnosis 
and prognosis of atherosclerotic evolution. 

1 0 Furthermore, regarding the above experimental results of the 

application of the EVE method to characterize human artery, while an 
optimal range of intraluminal pressures seems to be indicated to improve 
plaque detectability, the results also showed two specific features. Firstly, 
comparing elastography with histology, the geometry of the plaque 

15 appears to be preserved in the LSME elastograms (see Figures 25C and 
25D). For instance, the maximum plaque thickness measured in the 
elastogram of Figure 25C is close to 360 |j.m; this estimation is strongly 
supported by histology measurement conducted by Brusseau etal. (2001 ), 
who found a maximum plaque thickness of approximating 350 j^m for this 

20 very same carotid artery segment. Secondly, regarding biomechanical 
properties, a strain ratio close to 3 could be observed between the 
atherosclerotic plaque and the healthy surrounding vascular tissue for all 
elastograms presented in Figures 25. Such information may provide 
interesting insights about plaque components; it thus may help predicting 

25 plaque rupture and also help in therapy planning. The possibility of 
assessing quantitative strain values with the LSME represents an 
advantage of the present EVE method over the prior art. 
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A major advantage of the present EVE method over 
correlation-based techniques [de Korte et al., (1997; 1998; 2000a;) 
Brusseau et aL, (2001); Talhami et al., (1994); Ryan and Foster, (1997); 
5 Shapo et al., (1 996a)] stems from the fact that it allows computing the full 
strain tensor. For instance, complex tissue deformations such as rotation, 
scaling and shear can appropriately be assessed, whereas they are 
known to set a potential limitation for correlation-based methods. 

10 Although the present invention has been described 

hereinabove by way of preferred embodiments thereof, it can be modified 
without departing from the spirit and nature of the subject invention, as 
defined in the appended claims. 
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